Emergence of Topological Fermi Liquid from a Strongly Correlated Bosonic System in 
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Recent experiments on quantum degenerate gases give an opportunity for simulating strongly- 
correlated electronic systems in optical lattices. It may shed light on some long-standing puzzles in 
condensed-matter physics, like the nature of high-temperature superconductivity in cuprates that 
had baffled people over two decades. It is believed that the two-dimensional fermionic Hubbard 
model, or t-J model, contains the key to this problem; but the difficulty of unveiling the mystery 
of a strongly-interacting fermionic system is also generally acknowledged. Here, as a substitute, 
we systematically analyze the property of bosonic t-J model simulated in optical superlattices near 
unit-filling. In particular, we show the emergence of a strange topological Fermi liquid with Fermi 
surfaces from a purely bosonic system. We also discuss the possibility of observing these phenomena 
in ultracold atom experiments. The result may provide some crucial insights into the origin of high- 
Tc superconductivity. 

PACS numbers: 03.75.Hh, 03.75.Lm, 64.60.Cn, 74.20.-z, 74.20.Mn 



CD 



I 

O 

o 



> 

(N 
O 

o 



X 



In recent years, the physics community has witnessed 
a series of exciting discoveries and achievements. Among 
them, using ultracold atoms that form Bose-Einstein 
Condensates (BEC) or Fermi degenerate gases to make 
precise measurements and simulations of quantum many- 
body systems, is quite irnpressive and has become a 
rapidly-developing field [l|, Since atoms are cooled 
down to temperature near absolute zero and trapped in 
optical lattices building from six orthogonal laser beams, 
they provide us a peaceful playground for manipulating 
atoms with unprecedented accuracy. Some pioneering 
works Q revealed the promising potential of applying ul- 
tracold atoms to make quantum computer and quantum 
simulator: By changing the intensity, phase and polariza- 
tion of incident laser beams, one can tune the Hamilto- 
nian parameters including the dimension, the hopping 
strength and the particle interaction at will. People 
have successfully observed the Mott insulator-superfluid 
transition in both bosonic'?! and fermionic 'B!, Q degen- 
erate gases, and have demonstrated how to producep?] 
and control effective spin interactions in a double-well 
ensemble 0. All these evidences imply that an era in 
which atomic and optical physics unites with condensed- 
matter physics is within sight. 

Particularly, the two-dimensional fermionic Hubbard 
model (or t-J model) is one of the most interest- 
ing issues depicting the nature of high-temperature 
superconductivity d, [l^, [HI. Ever since the discovery 
of high Tc cuprates, tremendous efforts had been con- 
tributed to investigations of this model. Derived from 
Hubbard model at half-filling (one electron per site) , the 
t-J model describes the motion of doped holes in an an- 
tiferromagnetic (AF) background. It carries the essence 
of a strongly correlated electronic system with intrinsic 
competition between superexchange interaction of spins 



and hopping processes of charge-carriers (holes). Over 
the past two decades, people have employed many meth- 
ods and developed different schemes, hoping to fully un- 
derstand this model. Although a lot of consensus have 
been accomplished, there are still some ambiguities to 
be clarified. For example, can superconductivity evolve 
from a purely fermionic repulsive many-body system? 
Are there exist some new phases of matter in high-Tj, 
superconductors? 

Therefore, some people sug gest to use its bosonic coun- 
terpart, the bosonic t-J model]!^. [isl [T^ . as a trial model 
to investigate, for bosons are much more easier to deal 
with in both analytic and numerical approaches. After 
obtaining some experiences and conclusions, we may use 
them as reminders and analogies to original fermionic 
model. Besides, considering the present situation of ex- 
periments in ultracold atoms, where bosons are more 
accessible to be cooled and controlled, we believe it is 
worthwhile to explore the bosonic t-J model in optical 
lattices which can be formally written into two parts: 



H = Ht+Hj 



H.c 



(1) 
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where Hilbert space is restricted by the no-double- 
occupancy constraint J^a^la'^i'^ — ^' '^i<^ annihilates 
a two-component boson (cr =t,i denoting two internal 
states) and S.; is the (pseudo)spin operator at site i. Si = 
\ I]q/3 ^l^f^afiaip with Pauh matrix cr = {(7^^,(7^, cr^), {ij) 
denotes the nearest-neighboring counting. Ht and Hj 
describes hopping {t > 0) and an effective AF superex- 
change ( J > 0) interactions respectively. 

In a seminal experiment Bloch et al. loaded 
a two-component ^^Rb condensate, |_F = 1,tof = 1) (T) 
and \F — 1, mp — —1) (i), into arrays of isolated double 
wells. By changing the relative phase of the two laser 
standing waves, the potential difference A in one double- 
well can be raised or ramped down. The oscillation of 
the condensate after such manipulation contained infor- 
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FIG. 1: Illustrations of optical superlattices. (a) Construction 
of an optical superlattice. Combining two sets of optical lat- 
tices (wavelength A and 2A) with a relative phase difference 
TT to produce a superlattice. V is the depth of the poten- 
tial well. The superexchange processes are realized through 
virtual hoppings, resulting perturbations (7± A in the Hamil- 
tonian. (b) A visualization of a 2D superlattice, blue atoms 
denote one spin species (|T)) and red atoms denote the other 
species (|J,)). 



mation of the strength and sign of the effective superex- 
change energy J. With proper bias A, they could convert 
ferromagnetic interaction (J < 0) into antiferromagnetic 
one (J > 0). 

To implement a controllable hopping term, we have to 
generalize the double-wells to a set of biased optical su- 
perlattices. As illustrated in Fig. (IT|), for instance, we can 
start from the situation with one boson per site, differ- 
ent species in neighboring sites. Then by adding another 
set of optical lattice adiabatically, whose wavelength is 
as twice as the original one and has a relative phase dif- 
ference of TT, we can build a superlattice (with shallow 
sites and deep sites) , where a certain type of bosons be- 
ing trapped in a certain type of sites, for instance, ||) in 
shallow sites and ||) in its neighboring deep sites. 

By changing the shape of the superlattice, one can lift 
the potential difference A to a fixed value to realize the 
effective AF superexchange interaction (J > 0). Fur- 
thermore, by adjusting the relative phase and polariza- 
tion between incident lasers, one can control the tunnel- 
ing amplitude ti^ of different inner states of atoms in 



different sites to be identical. Therefore, we can intro- 
duce the hopping effect while keeping a(n) (inhomoge- 
neous) Heisenberg-type interaction in an optical super- 
lattice, (see Methods for details.) 

In addition, according to ref. [l^, one can introduce 
vacancies in a BEC trapped in one-dimensional (ID) op- 
tical lattice by pointing an electron beam at specific sites 
to remove atoms. Due to the small diameter (100 ~ 150 
nm) of the electron beam compared with the lattice spac- 
ing, holes can be doped in without additional losses in 
other sites. The final occupying configuration can be vi- 
sualized with highly spatial resolution through scanning 
electron microscopy. Since single-occupancy has been 
achieved in current experiments [sl. [l^ as an initial setup, 
we believe that this technique can be further extended to 
2D and 3D superlattices loaded with multi-component 
atoms, thus one can manipulate and(or) remove different 
species deliberately and observe the dynamics of holes 
with single-atom, single-site sensitivity. Above all, the 
realization of theoretical model ^ as well as the nec- 
essary doping technique are all within reach in current 
experiments. 

The main analysis is arranged as follows: for bosonic 
t-J Hamiltonian near unit-filling, we analyze the doping 
effect in a topological perspective. We show the existence 
of several exotic phases as vacancies being gradually in- 
troduced. Then we discuss how to detect these features 
in currently-available techniques in ultracold atom exper- 
iments. 



I. HOLONS: BOSONS OR FERMIONS? 

At unit-filling, the bosonic t-J model is simple and can 
be reduced to AF Heisenberg model, H = JJ2(ij) ' ^j- 
A variational wave-function based on Resonating Va- 
lence Bond (RVB) picture by using Schwinger-boson de- 
scription can produce [13, E, 13 ^'^ unrivaled accurate 
ground-state energy; a generalized versionpjj] can further 
provide the staggered magnetization and spin excitation 
spectrum precisely. So the bosonic RVB picture is a nat- 
ural choice for describing AF Heisenberg model. Explic- 
itly, we can introduced two flavors of bosons on each site, 
created by a canonical operator foj^ acting on the vacuum 
|0) without bosons, satisfying bl^bia = 1. This represen- 
tation is equivalent to the following operator identity be- 
tween the spin and boson operators Si = ib^abi. Here 

bi = (6i|, fcix)"^ is ^ bosonic spinon annihilation operator. 
The mean-field value is characterized by a bosonic RVB 
order parameter Af^ — {bicrbj,-cr) 7^ for the nearest- 
neighbor sites, which depicts the short-range AF corre- 
lation as (Si ■ Sj) = — i|Afjp. At zero temperature, 

spinon b becomes massless and Bose-condensation takes 
place with (b) ^ 0, corresponding to the long-range Neel 
order in x-y plane. 

To learn the property of an AF order with vacan- 
cies, we should generalize the idea of spin-charge sep- 
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aration, which has become a very basic concept in un- 
derstandin g th e doped Mott insulator related to high-Tj, 
cupratefiol l20l|. Unlike a usual quasi-particle that carries 
both spin and charge quantum numbers in conventional 
metals, it states that the system has two independent el- 
ementary excitations, the neutral spinon and the spinless 
holon. It is assumed that usual quasi-particle excitations 
may no longer be stable against the spin-charge separa- 
tion mechanism once being created, e.g., by injecting a 
bare vacancy into the system; and it has to decay into 
more elementary spinous and holons. In other words, 
to theoretically describe the introduction of a single va- 
cancy, one has to first annihilate a particle state j^") to- 
gether with a bosonic spinon bia-, then generate a spinless 
operator hi denotes a holon (a vacancy) as 



(2) 



In addition, hi and bia should satisfy the no-double- 
occupancy constraint, hi hi + J2a ^la^iiy = 1- In contract 
to the case of high- Tc cuprates, holons in bosonic t-J 
model are neutral for they are actually the absence of 
cold atoms in a certain site. 

However, the situation changes when the hole moves. 
According to Marshall 2l| , the ground-state wave func- 
tion of the Heisenberg Hamiltonian for a bipartite lat- 
tice satisfies a sign rule. It requires that flips of two 
antiparallel spins at nearest-neighbor sites are always 
accompanied by a sign change in the wave function: 
I- • • ti • • • ) '-^ (-1) X I- • • it • • • )• To show the sign effect 
in detail, we divide a bipartite lattice into odd {A) and 
even {B) sublattices and assign an extra sign (—1) to ev- 
ery down spin at A site. When the hole initially locating 
at site i hops onto a nearest-neighbor site j, the Marshall 
sign rule is violated, resulting in a string of mismatched 
signs on the vacancy's course[22i . The spin wave function 
is changed into 



I*) 



where (—1)^ = 1, (—1)''' = —1- The phase-shift factor 
e'® can also be regarded as a unitary transformation on 
an arbitrary operator: O s- e*®Oe~'®. For example, op- 
erators of holons hi and spinous bi^ are transformed as 
follow. 

Now we have restored the simple spin wave function j'l') 
as equation ([2|) originally defined, but have got a set of 
nontrivial operators, and it is just the aim of the above 
transformation. ^ 

Furthermore, by defining h'^ = /jte"' , we can 

see that except for a phase factor , holon hi and 

spinon bi^ are symmetric. 



(3) 



h — >h, i- 



h — >h. ia- 



namely, hi ^^ — > bi^ and 6 jo- ^ — ^ hi. The phase 
factor e 2 ^ j means an additional lattice 7r-flux-per- 
plaquette for holons. From equation ([3]), one can see 
that there exists a mutual semionic statistics between 
holons and spinous, where holons perceive spinous as tt- 
vortices and vice versa. In particular, after the trans- 
formation h'^ = /i|e~* ^'^"j , h'^ obeys fermionic anti- 
commutation [2^ , 

{h'm}^5,,. 

Finally, based on the RVB ground state by considering 
the phase string effect [23, 24, 25], the effective Hamilto- 
nian of holons in bosonic t-J model can be written as 



A 



■'oh'^^h'^+H.c. 



(4) 



where njg^ | — ^Jg^ j&iGA.i is the number of down spin 
on each yl-sublattice. In particular, if the hole moves 
through a closed path C on the lattice to return to its 
original position, it will get a Berry phase (— 1)^^^ where 
is the total number of down spins "encountered" by 
the hole on the closed path C. [1^ [H, [IH] This process is 
illustrated in Fig. ((2^). 

To deal with this unavoidable Berry phase l)^c 
in the ground-state wave function when there is mo- 
bile hole, we introduce a phase-string transformation[23j 



I*) 



e*e|*) = I*), where Q ^ 



n 



3l 



and 

are occupation number operators of the hole and 
down-spins, with a phase factor =Im[ln(zi — Zj)]. Here 
z = X + iy denotes position and the subscripts i and j 
denote lattice sites. Considering the single-occupancy 



constraint, 9 = — ^Ey^i'^ufl 



-1 



n 



where ~ l^fjl* is the effective hopping amphtude 
of holons, fj, is the chemical potential. The gauge 
field Qij satisfy the topological constraint J^c '^v ~ 
±7r^;gp(nf| — n^) for a closed loop C, 4>^j describes 
a TT flux per plaquette, 0? — ±tt. It is obvious that 
the difference between fermionic t-J model and bosonic 
t-J model is the statistic of holons: in former case, holons 
are bosons; and in latter case, they are fermions. 

Now we may be able to answer the question: what is 
the fate of a holon in an AF order? As we have mentioned 
above, the AF order lying in x-y plane can be seen as 
a Bose condensation of spinous, {bia) i^ 0. According 
to the above effective Hamiltonian and the quantity 
Sj = (— l)*5j|6j|e*^j , a holon introduced here is a 
topological defect carrying spin twists and changing its 
peripheral spin configuration, as a meron-like obiect[27l. 
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m, [2^ [sOl , such a spin vortex (for a holon at the origin) 
may be characterized by a unit vector = (—1)* (S^) / S, 

This meron-hke spin configuration is schematically shown 
in Fig. (H)). 

By now, one may imagine that a single hole would 
become a meron in a long-range AF order. However, 
the answer is not quite right. For a single meron 
configuration, the energy is logarithmically divergent, 
E w J|A|^ |ln (L/ao) with L the size of the system. 
In order to remove this infinite-energy cost, as in the 
present case, each holon-meron has to "nucleate" an 
anti-meron from the background spontaneously. De- 
fine nf -I- in^ = e"^""'"*"^' , with the unit vector ng = 
(cos 00, sin</)o) as the magnetization direction at infin- 
ity. In presence of a fermionic holon centered at —Zk/2 
and an anti-meron centered at we havefsTi [3^ 

= Imln [{zi - Zk/2) / {zi + Zk/2)] with Zk = el + ie^. 
By using Gfc = (e|, e^) to denote the spatial displacement 
of the holon and anti-meron centered, a dipolar spin con- 
figuration at a sufficiently large distance is obtained, 

« --2 ' > l^fcl- 

Thus, each pair of holons and anti-merons forms a com- 
posite, as shown in Fig. ((^t). In contrast to the loga- 
rithmically diver gent meron-energy, the energy a dipole 
Ed becomes finitelHlll, Ea ~ J|Afj | In [(|efc| + «„) /ao], 
|efc| > ao- 

Physically, one may consider a bare hole (spinless 
holon) created by an annihilation operator at point 
a, then it may jump to point b via some discrete steps, 
being connected by a phase string in between. Since the 
holon can reach h through different virtual paths orig- 
inated at a, this singular phase string is then replaced 
by, or relaxes to, a smooth dipole configuration. An anti- 
meron itself is a semi- vortex formed by condensed spinous 
and it is immobile. Therefore, the hole-dipole as a whole 
must remain localized (self-trapped) in space. The re- 
sulting spin configuration has a dipolar symmetry while 
the distortion of the direction of magnetization is long 
ranged and decays as r~^[3^[33t. 

II. TOPOLOGICAL FERMI LIQUID 

In the lightly doped region, (5-^0, there exist local- 
ized holes that are self-trapped around anti-merons via a 
logarithmically confining potential V {r) k, g^ln(|r|/a), 

= 27rJ|A^j|, and their dipolar moments are randomly 
distributed. Such kind of AF ordered state with random 
dipoles has been studied in refs. [Sll, [H, Is^, HI]. It is 
known that when T = 0, the spin-correlation length is fi- 
nite for arbitrary doping 5. This implies the destruction 
of long-range AF order, as long as (5 > 0. Consequently, 
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FIG. 2: (a) A schematic demonstration of the phase string 
effect on a long-range AF order background. The point in 
A (i3)-sublattice is denoted by red(blue) arrows. The double- 
dot path is the course that a hole (black open circle) moves 
through. For clarity, spins on this path is plotted in light 
colors. Once the hole's hopping results a down spin on A- 
sublattice, there is a sign change ( — 1) over the wave function 
j^), as described in the text, (b) An illustration of a meron 
configuration. Due to the mutual statistics, the vacancy has 
become a fermionic particle, (c) An illustration of a holon- 
dipole configuration. It is a confined object composed of a 
holon-meron and an anti-meron at two poles: a(— 2^/2) and 
b{zk/2), connected by a branch-cut with the spatial separa- 
tion Gfc as a dipole moment. 
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AF order is limited mainly in finite sizes, where the size 
^ of a domain is determined by the hole's concentration, 
^ w a /Vs. Accordingly, this state has been termed as a 
cluster spin glass (39l[40|. The spin glass freezing temper- 
ature is then expected to vary as Tg a f ^ a 0^/(5, below 
which holons tend to form a glass and their dynamics 
strongly slows down. 

When doping increases, more and more dipoles ap- 
pear. As a result, the "confining" potential V (r) be- 
tween these dipoles will be screened by polarizations 
of pairs lying between them. When this screening ef- 
fect becomes so strong that even the largest pair has to 
break up, holons and anti-merons are liberated to move 
individually. This qualitative analysis suggests that a 
localization-delocalization phase transition should occur 
at a critical point Sc, in a fashion of the Kosterlitz- 
Thouless (KT) transition. Using the renormalization 
group method [sH-Iss!]. one can determine the critical hole 
density 5c = Sc(T) or temperature Tde = Tdc((5) at which 
dipoles collapse and holon-merons are "deconfined" from 
the bound state with anti-merons. At zero temperature, 
the critical hole density has been numerically determined 
as 

6ciT = 0) « ^ = 0.043. 

Therefore, a zero-temperature quantum critical point ex- 
ists at Sc — 0.043 where hole-dipoles dissolve into holon- 
merons and anti-merons. 

When S > 0.043, fermionic holons are delocalized, so 
a Fermi surface emerges in the momentum space. And 
since two holons always repulse each other, for they 
are merons with same topological charges, then these 
fermions cannot pair with each other and it makes the 
Fermi surface stable. This is different from the case where 
electrons form Cooper pairs in superconductors, because 
here we have topological repulsions not electronic repul- 
sions. 

Thus, we display the global phase-diagram of bosonic 
t- J model in Fig. (3). There exists three different regions 
globally: region I, AF insulator, 5 = 0; region II, insu- 
lating spin-glass (SG), S e (0,0.043); region III, topo- 
logical Fermi-liquid (TFL) phase with spin-liquid ground 
state, S G (0.043,0.5). The phase boundary that pro- 
tects the physics in low doping regimes can be written 
as fcsT « c(l — 2S) J, c is an integration quantity (see 
Methods). Particularly, our results illustrate the emer- 
gence of an exotic fermionic state (TFL) in a pure bosonic 
system. 

In TFL state, the effective Hamiltonian of holons can 
be written as 

(ij) i i 

and the dispersion relation is 

Ek=^i± [cos2(fc,) + cos2(fcj,)] + A2, 




FIG. 3: The global phase diagram. The solid black line indi- 
cates the phase boundary where |A|j| = 0. Different regions 
are demonstrated in the text. The inset depicts detailed phase 
boundaries as well as the quantum critical point 5c in low dop- 
ing region (5 < 0.1). The Fermi energy curve ^ {5) implies the 
existence of Fermi surface of deconfined holons. It is scaled 
by the superexchange energy J as shown in the inset, where 
we have chosen U/t = 12, namely, tj J = ?> and A/J = 50. 

as Fig. shows. 

In Fig. ([5]), we show the evolution of Fermi surfaces 
by plotting four doping cases, 5 = 0.05, 0.1, 0.15 and 0.2, 
respectively. From the global phase diagram Fig. (3) 
and Fig. ([5]) , we suggest that in experiments people may 
observe TFL in bosonic t-J model of 15% hole concen- 
tration, of which highest Fermi energy with clear Fermi 
surfaces are predicted. 

Besides, when holon moves, the spin configuration 
changes simultaneously, thus there is no long-range AF 
order anymore. Instead, the magnetic ground-state be- 
comes a spin liquid state with invariant spin-rotation and 
translation symmetry. Without anti-merons matching 
with holons, holons exert an "effective magnetic field" on 
spinous (recall the mutual relation between the two exci- 
tations), Bh = T^Ph with ph the density of holons. Hence, 
a new length scale is introduced to the spinon system, 

— 1/2 

which is the magnetic cyclotron length Ic = Bf^ ' . Ic 
will later be connected to the remaining magnetic corre- 
lation length. 

As we mentioned above, doping creates a Landau-level 
structure in the spinon spectrum. Thus, low-lying spin 
fluctuations are expected to be sensitive to dopants. Af- 
ter calculations, we find 

X" (q,c.)^x" (Qo,c.)e-l'i-Q«l^'>. 
Here x" (Qj^) is defined as (/3 = l/fc^T) 

X" (q, ^) = ^ (1 - e-^^) J dtdve^^'^'-^-^^ (S (r, t) ■ S (0, 0)) 

(5) 
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FIG. 4: The dispersion relation of the topological Fermi liq- 
uid. The relative energy Ek — fJ, of the upper branch is shown 
in the unit of t with A/t = 20. 
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FIG. 5: Illustrations of Formi surfaces. The doping concen- 
trations are shown as titles of individual plots. 



where (S (r, t) ■ S (0, 0)) is the spin-spin correlation func- 
tion. This is a Gaussian type around the AF wave- 
vector Qo = (±7r, ±7r)/ao with a width l^^. Conse- 
quently, it determines the spin-spin correlation in the 
real space as cos (Qo • r) e"'""' , where the correlation 
length ^ ~ V^lc — ao-v/2/ (ttS) is in the same order of 
the average hole- hole distance[41j. Namely, doped holes 
break up the long-range AF correlation into short-range 



AF fragments within a length scale of ^. 

In contrast, the ground-state of fermionic t-J model 
is a superconducting state of condensed holons, in 
which spin liquid state is always masked by holons' 
condensation ^42.] . While in bosonic t-J model, it is the 
topological Fermi liquid state with massless fermionic ex- 
citations, in which one can easily detect properties of this 
spin liquid state. For these reasons, the observation of 
this exotic quantum state (TFL with short range spin 
correlation) in bosonic t-J model may pave an alterna- 
tive approach to verify the microscopic theory of high Tc 
superconductivity. 



III. DETECTIONS AND SUMMARY 

In region I and II, the system is an insulator in the 
first place, so the momentum distribution is universal and 
isotropic in the Brillouin zone (without considering the 
trapping and boundary conditions). This feature can be 
detected by measuring density profile measurement [iij , 
or via Bragg spectroscopy [43] to observe the momentum 
spectrum. Additionally, in the Neel state (I), the con- 
densate exhibits long range AF order, while in the spin 
glass state (II), this long ranged magnetic order disap- 
pears. This difference can be distinguished through spa- 
tial spin-spin correlations, {Sz (ri) Sz {r2)) ■ As demon- 
strated in refs.ji^ H^, using a probe laser beam which 
goes through the condensate, one can measure the phase 
shift or change of polarization of the outgoing beam 
(A'out) to obtain the magnetization (Xout) oc (Mz) cx 
J drcf) (r) {Sz (r)), where (r) is the spatial intensity pro- 
file of the laser beam. The quantum noise {X^^^) cx 
(Af2) cx / dridr^^ (n) (ra) {Sz {ri)){Sz (ra)) reveals 
the atomic correlations in the system. This method can 
directly examine the existence of 2D AF correlations. 

In region III, we have demonstrated that one natural 
characteristic of TFL is the existence of Fermi surfaces 
in the Brillouin zone. Esslinger et al.fi^l had successfully 
observed the Fermi surface in a 3D optical lattice filled 
with fermionic atoms. Similarly, one may also observe 
these Fermi levels in our bosonic system, as shown in Fig. 
([5]). Besides, the short-range AF fluctuations may also 
be observed. As equation ([5]) indicates, the spin dynamic 
structure factor as well as the dynamic spin susceptibility 
function arc linked with spin-spin correlations. The latter 
one reflects the effective short-range magnetic correlation 
length ^. 

To sum up, in this paper, we first demonstrate the 
implementation of bosonic t-J model in optical super- 
lattices filled with a two-component BEC. And for the 
first time, we systematically discussed the possible quan- 
tum phases of this model upon doping. A key feature 
is the mutual semionic statistics between the two ele- 
mentary excitations of the system, holons and spinous. 
When there are only a few holes around, they tend to 
form hole-dipoles, locally changing the underlying spin 
texture into meron-anti-meron pairs. At low doping, the 
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deformed spin configuration loses the long-range order, 
and there emerges a spin glass state. When holes are pre- 
vailing, through a quantum phase transition {6c — 0.043) 
that frees those self-localized holons in spin glass state, a 
strange topological Fermi liquid with Fermi surface ap- 
pears in a purely bosonic system, as a significant result 
of the mutual statistics between holons and spinous. We 
also list accessible experimental approaches to verify and 
detect these novel theoretical predictions. 



detailed terms, we can obtain the effective Hamiltonian 

as 

^cff = -t^{al^aja + H.c.) + 

Y,[J^s^s]~J±is-s^ + s!sy)], (7) 

where 



IV. METHODS 

A. Perturbation Theory in Large U Limit 

The two-component Bose-Hubbard model in a biased 
superlattice can be generally written as 

H = -tY^ [al^a-ia + B.C. + ^ {h,a ~ nj^)] + 



U' 



^iT"-a + E "'"^ ^ ' 



(6) 



where t = ^le/irEr {V/Erf^'^ e^V^/^- j] with atomic 
recoil energy E^. — h?k'^/2m, k is the wave vector 
of laser, m is the atomic mass, V describes exter- 
nal potential. The inter-species repulsion is defined as 
U = yJ'&l'nkasEr iy jErf'^'^ with as the s-wave scatter- 
ing length among species, and the intra-species repulsion 
is V — ^JWpKka'gEr{V/Er)^^'^, a'^ is the correspond- 
ing scattering length. This term vanishes at exact unit- 
filling, but can contribute in higher-order tunneling pro- 
cesses. 

In the large U limit {U ^ t) and nearly unit-filling 
n < 1, which is actually a prerequisite for building the 
aforementioned superlattices and is naturally satisfied, 
the hopping term can be treated as a perturbation. Thus, 
we introduce two projection operator: V projects the 
initial Hilbert space onto the single-occupancy subspace 
(|0),;,|T),,|i),), 2 = 1-7' projects onto the double- 
occupancy subspace (|Ti)i)- We further divide equation 
dS]) into Hq + Tmix, where Hq describes processes in V 
subspace and Tmix mixes the upper and lower band via 
virtual tunnelings. 

Applying canonical transformation Hctf — e^^He^ 
and eliminating first order terms of T,„ix, we find 



^^cff ^ Hf) + - 




n) {n\ 



The projector |e) (e| can be chosen as V and Q. Energy 
differences between upper and lower band are Eq — E-p w 



?7±A,^±Athus5= {if- 
U = U~A^/U, il' = U' - 



^+2U'-^){VTQ-QTV) with 
A^/C/'. After calculating the 



4*2 
IT 



8^ 



4*2 
If' 



When as and a'^ do not differ very much (this condition 
can somehow be realized in experiments) , 



__4*2 

u 



with e denoting a small spatial inhomogeneity. In fact, 
this small inhomogeneity is necessary in our following dis- 
cussions, for it can be seen as an analogy for the effective 
inter-layer-coupling in high-Tc cuprates. 

At the first glance, it seems that bosons (say, |t)) in 
shallow sites may hop to neighboring deep sites more 
easily due to lower energy barrier they feel, comparing 
their counterparts (|i)) in the reverse process. However, 
thanks to spin-dependent controlling techniques, which 
can be achieved by adjusting parameters of incident light 
beams, one can regulate tunneling amplitudes of differ- 
ent species, making atoms in different sites feel a similar 



hopping matrix element, t| 



t. This is the reason 



we straightforwardly set an identical hopping term t at 
the beginning. 

Meanwhile, to achieve an effective AF superexchange 
interaction, we need to set A = ^/2U {U — —U) to 
change signs of coefficients: Jz ~ 'ti'^/U+e, J± = —At^/U 
in equation ([7]). As a result, the global coefhcient of spin 
interaction becomes positive, J ~ At'^/U, as equation ([T]) 
asks. 

Under these circumstances, we finally derive the 
bosonic t-J model with a small inhomogeneity in z di- 
rection, which can be formally written as equation (fl]l. 
Quite recently, two papers [i^, suggested a similar 
proposal and made detailed discussions. 



B. Mean field calculation from slave-particle 
approach 

This representation is equivalent to the following op- 
erator identity between the spin and boson operators. 

4* ■(^string ^ " t'-o.^ '■j. 

Defining aia- — hle^^" bia and Si — ^b-crbi, where hi 
creates a fermionic holon, leaving a non-local phase string 
and foio- annihilates a bosonic (Schwinger-boson) spinon 
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at site i, we rewrite equation ((T|) as 7? = Hh + Hg'. 

{ij)aa' 

Here 6"^""^ = - a^'^) is a topological phase 

with the contribution from spinon number fi^, <i&^ = 
Zli^^j.o- ~ ^Ol^fcr and the holon's {n'^) 
contribution, = [^^ i^i ~ ^i)] ■ ^i^'i ^{j a-ud 

A^j describe quantized fluxes bounded to spinous and 
holons. We then introduce a mean-field bosonic RVB or- 
der parameter A'' — ^„{e~^'^^^ihiahj^-a-), and take the 
constraint on total spinon number and dilute effect into 
account, we apply mean-field approximation on H and 
solve two self-consistent equations: 

5 = 2-l^Afci?fc-icoth(/3i?fc/2), 
fc 

k 



to obtain the phase boundary 



keT = 2J{1~25) [2-5) {\n\ 



that separates A" = and A** ^ 0. Here N is the total 
lattice number. Detailed definitions can be found in ref. 
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